Interaction of Mesoscopic Magnetic Textures 
with Superconductors. 
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We analyze magnetic textures in a thin magnetic film and screening currents induced by these 
textures in a superconducting film. The two films are supposed to be parallel and close each to 
other, but interact only via magnetic fields. We consider also vortices inside superconducting film 
and their interaction with magnetic textures. As examples of magnetic textures we consider magnetic 
dot and magnetic vortex. We derive the condition at which spontaneous formation of a vortex or a 
coupled vortex-magnetic defect is energetically favorable. It is proven that the normal component 
of magnetic field generated by a magnetic dot changes sign in the presence of the superconducting 
vortex. 
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The fabrication and experimental study |l]] of meso- 
scopic heterogeneous magnetic/superconducting systems 
together with recent theoretical predictions open a 
new class of physical effects. Earlier two of us (I.F.L. and 
V.L.P.) proposed to separate superconductivity and mag- 
netism in space employing the modern technique of nano- 
fabrication Qj. The proximity effect which suppresses 
both order parameters can be easily avoided by growing 
insulator oxide layers between ferromagnetic (FM) and 
superconducting (SC) components. Inhomogeneous mag- 
netization of the magnetic texture generates magnetic 
field penetrating into the superconductor. The magnetic 
field from the SC currents interacts with magnetic sub- 
system pM. We have proposed different realizations of 
mesoscopic magneto-superconducting systems: arrays of 
magnetic dots on the top of a superconducting film j|, 
magnetic/superconducting bi-layers 0, magnetic nano- 
rods embedded into a superconductor B. 

In the majority of proposed systems a magnetic texture 
interacts with the superconducting current. An inhomo- 
geneous magnetization generates magnetic field outside 
the magnet. The magnetic field generates screening cur- 
rents in superconductors which, in turn, change the mag- 
netic field. The problem must be solved self-consistently 
. In this article we develop a general formalism for inter- 
acting inhomogeneous magnetization (texture) and su- 
perconductors in the London's approximation. Employ- 
ing this formalism, we find two elementary solutions for 
a circular magnetic dot on the top of a superconducting 
film, magnetized perpendicularly to it, and for a cou- 
pled magnetic and superconducting vortices. London's 
approximation works satisfactory since the sizes of all 
structures in the problem exceeds remarkably the coher- 
ence length £. 

Recent theoretical works consider problems close 
to ours. The main difference between our work and the 



Ref. H is that they do not incorporate spontaneous vor- 
tices in the ground state. In comparison to the Refs. 

our approach is more general and universal. A more 
detailed comparison with these references will be given 
later. 

So far only arrays of sub- micron size magnetic dots cov- 
ered by thin superconducting films have been prepared 
and studied experimentally 0]**. The effect of commen- 
curacy on the transport properties was reported. This 
effect is not specific for magnets interacting with super- 
conductors and was first found many years ago by Mar- 
tinoli and his group II. Predicted effects specific for 
FM/SC systems, including topological instability of ho- 
mogeneous state in bi-layers 0j, were not yet observed 
exp er iment ally. 

General formalism The total energy of a stationary 
magnet-superconducting system reads: 
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where B is the magnetic induction, M is the magnetiza- 
tion, n s is the density of superconducting carriers, v s is 
their velocity and m is their effective mass. We assume 
that the superconducting density n s and the magnetiza- 
tion M are separated in space. We assume also that the 
magnetic field B and its vector-potential A asymptoti- 
cally turn to zero at infinity. Employing Maxwell equa- 
tion VxB = — j , and B = V x A, the magnetic field 
energy can be transformed as follows: 
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The current j can be represented as a sum: j = j s + j m of 
the superconducting and magnetic currents, respectively: 
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In the gauge invariant Eq.|| ip is the phase of the SC carri- 
ers wave-function and <fio is the (SC) flux quantum. Plug- 
ging Eq.([|) into Eq.(||) and Eq.(|l|), after minor trans- 
formations, we arrive at following equation for the total 
energy: 
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where a = A(r, z = 0) and b = B(r, z = 0). The vector- 
potential satisfies Maxwell-London's equation: 

V x (V x A) = — -AS(z) + —VipS(z) (6) 

A mc 

+ 47rV x (mS(z)) 

where A = A|/ds, is the effective screening length for the 
SC film, Al is the London penetration depth and d s is 
the SC film thickness. 

According to our arguments, the term proportional to 
Vip in Eq. (0) describes vortices. A plane vortex char- 
acterized by its vorticity q and by position of its center 
on the plane ro contributes a singular term q z *Jf~*^ to 
V<^> and generates a standard vortex vector-potential: 
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Different vortices contribute independently into the 
vector-potential and magnetic field. Below we calculate 
the contribution of a magnetic texture to the vector po- 
tential A m (R) in a special gauge A z = 0. The magneti- 
zation vector field can be represented as follows: 



m(r) = m z (r)z + m"(r) + m ± (r) 



(8) 



where m"(r) is the gradient part of m(r) and m J -(r) can 
be represented as a vector product z x V/ ( / is an arbi- 
trary function of r ). This representation is unique. The 
corresponding contributions to A m are : 



A m (R) = A"(r) + A ± (R) 



(9) 



where 



The phase gradient can be included into A as a 
gauge transformation with the exception of vortex lines, 
where tp is singular. Eq. (Q) allows to separate the en- 
ergy of vortices from the energy of magnetization induced 
currents and fields and from the interaction energy. 
Two dimensional textures and vortices. Below we apply 
our general formalism to the case of two parallel films, 
one FM, another SC, both very thin and very close to 
each other. Neglecting their thickness, we assume that 
both films are located at z = 0. The super-carrier den- 

( 2 ) 

sity n s (r) can be represented as n s (R) = 5{z)ni>(r) and 
the magnetization M(R) can be represented as M(R) = Here : 

(2) 

<5(z)m(r), where nl (r) is the super-carrier density per 
unit area, m(r) is the magnetization per unit area and 
r is the two-dimensional radius-vector. In what follows 
is assumed to be a constant and the index (2) is 
omitted. The energy (4) can be rewritten for this special 
case: 
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where non-zero component Gij(R — r') are (R = (r, z)) : 
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e a g is the second rank antisymmetric tensor and J 7 (x) 
are the Bessel functions. 

Equations (jl^,[l2|) can be further simplified if the mag- 
netization has only z and r components and both depend 
only on r. Then the vector-potential A™(R) is directed 
along azimuthal lines and its only non-zero component 
reads : 



A£(R) =J2 [ K a (R,r')m a (r')d 2 r' 



where a = stays for z, and a — 1 stays for r, 
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and 
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Santos et al. ^ have developed a formalism for calcu- 
lation of magnetic fields and screening currents generated 
by magnetic textures in superconductors similar to ours. 
However, they did not consider singular current distri- 
butions, i.e. vortices. Bulaevsky et al. considered a 
special system of prepared domains in a thick magnetic 
slab contacting with a semi-infinite superconductor and 
generating vortices in it. The spontaneous formation of 
domains and chains of vortices was discussed earlier [|To"P] 
for a bi-layer composed from very thin films. 

Magnetic Dot. For infinitely thin circular magnetic dot 
with magnetization m = mz and radius R the magnetic 
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field can be calculated using Eq.|L3|Jl 
has the form: 



Its z-component 
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Numerical calculation based on Eqjlf show that B z on 
the film (z — 0) changes sign at some r > R if q is pos- 
itive (see Fig.Q), but it is negative everywhere at r > R 
if q = 0. The physical reason for this behavior is simple: 
the dot itself, without vortex currents, generates dipolar 
field which has the sign opposite to dipolar moment in 
the plane passing through the dipole and perpendicular 
to it. The SC current resists this tendency. At small 
distances from the center the screening is neglegible and 
the field generated by dot dominates. At large distances 
r ^> A the field generated by a vortex decays like 1/r 3 , 
whereas the screened field of the dot decays as 1/r 4 . The 
position of node ro can be easily determined if R <C A: 
then ro = 2nyJ mi? 2 A/$o- Thus, the measurement of 
magnetic field near the film may serve as a diagnostic 
tool to detect SC vortex bound by the dot. The experi- 
mental measurements Jn| demonstrated that B z changes 
sign, i.e. the dot indeed generates the vortex. 

In the presence of a vortex with charge q , the energy 
of the system can be calculated using Eqs.(§, [if). The 
result is: 
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where e m = m<I>o is the characteristic SC/FM inter- 
action energy. The vortex becomes energy favorable if 
e m R/(2\) > e or mR > [$ /(8tt 2 )] ln(A/f). The ratio 
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In 



is the relative strength of the SC/FM coupling j 
Eq. @ n m is the density of magnetic atoms, d„ hS are 
the thicknesses of magnetic and superconduccting films, 
S is the value of an elementary spin in the magnet and g 
is the Lande factor. Far from the threshold the vorticity 
q generated by the magnetic dot is an integer closest to 
the value 
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The logarithmic factor in Eq. ([jl]) varies in the range 
3-8. For typical values R — 10 _4 cm, d m — 10 _6 cm, 
and n m — 10 22 cm~ 3 , we find q > 1. Eq. (19) shows 
that, apart the logarithmic factor, the SV generation is 
controlled by properties of the magnetic film. The char- 
acteristic scale in superconducting film is A. For this rea- 
son, at R ~ A, the condition S ~ 1 gives an approximate 
criterion for the spontaneous vortex formation. 
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FIG. 1. Normal to the dot magnetization as a function of 
distance from the dot's center measured in the units of dots 
radius R. Top: without vortex, for A = 20R; Bottom: with 
vortex the case when q0o/(47r 2 m) = 4 and A = 207?. 

Sasik and Hwa Q mimicked the magnetic dots by mag- 
netic dipoles. This approach does not reflect properly the 
real dot geometry. As it is seen from our results, the mag- 
netic field and energy can not be expressed in terms of 
the total dot magnetic moment (irmR 2 ) only. 
Coupled Magnetic/Superconducting Vortices. Earlier we 
predicted spontaneous formation of the pairs of bound 
superconducting and magnetic vortices in the case of bi- 
layer consisting from superconducting and magnetic film 
with easy plane anisotropy H . A detailed analysis given 
below confirms this prediction, but reveals more restric- 
tive criterion for this phenomenon. The exchange Hamil- 
tonian of the magnetic layer reads: 



H M = J -f J d 2 r(Vn) 2 



(20) 



where n(r) is the unit 2D vector directed along the lo- 
cal magnetization. Further we assume that the vector n 
lays in plane (XY symmetry.) In a single magnetic vor- 
tex (MV) the local magnetization is directed along the 
radius: n(r) = r/r. The energy (20) is invariant under 
a global rotation of all spins by the same angle. The 
MV energy grows logarithmically with the system size 
L : E w ttJm ln(L/£). We employ general formulae, (0, 
|IJ,|l5|) to calculate the vector-potential, magnetic field 
and the total energy. The vector-potential induced by 
the vortex, including the effect of screening currents is: 
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Employing the standard vortex vector potential Eq. |?] 
with R = and eq. (||), we calculate the energy of the 
MV-SV pair with logarithmic accuracy: 



E = q 2 e ln( 



^)-^ln(^) + 8^m 2 Aln(^) (22) 



The exchange energy of the magnetic vortex E ex is typi- 
cally by 1-2 order of magnitude smaller than eo and can 
be neglected. The energy is minimal at q being an in- 



teger part of the value q 
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The coupled pair 



of SC/FM vortices becomes energy favorable when mini- 
mum is realized at q = 1 . It means that S must overcome 
a threshold value S r — in ^^\ ■ In the limit L — > oo, 

ln(L/A) 

due to the logarithmically divergent second term in Eq. 
|22| , the magnetic/superconducting vortex pair is always 
energetically favorable. However, in practice the loga- 
rithmic factor never exceeds 5 — 7. Therefore, the value 
5 exceeding 1 triggers generation of the MV/SV pairs. 
The value of 5 can be enhanced by either reducing the 
supercarrier density n s or increasing the thickness of the 
magnetic film. 

In a bi-layer of XY ferromagnet and SC films with 
Curie temperature T c larger than SC transition temper- 
ature, the ferromagnetic state will be observed in the 
temperature interval between T s and T c . Below T s the 
MV/SV vortices proliferate. They destroy both, ferro- 
magnetism and superconductivity. Hovewer, the reen- 
trant phase transition at T = Ty leads to restoration of 
both order parameters if 8(T = 0) is smaller than the 
logarithms ratio. 

The vortices of the same sign repulse each other, 
whereas the antivortex attracts vortex. Together they 
form a dense system, most likewise a crystal at low tem- 
perature, with the lattice constant d ~ A, since A is the 
only characteristic length in the system. This is a new 
phase, never observed in 2d systems before. It differs 
from a vortex crystal in the 3d superconductors in mag- 
netic field. It is also very different from the Berezinskii- 
Kosterlitz-Thouless zero-field liquid vortex phases having 
no bound vortex-antivortex pairs. 

In conclusion, we presented a general formalism to 
treat the interaction between magnetic textures and su- 
perconductors. As applications, we have shown that vor- 
tices in superconducting films can be generated by mag- 
netic dots with normal to the film magnetization. Su- 



perconducting vortices also appear in bound pairs with 
magnetic vortices in a homogeneous, easy-plane magnetic 
film. We have found how the magnetic dot size and 
its material control the vortex appearance. Spontaneous 
generation of the MV/SV pairs in the bi-layer of the XY 
magnetic and SC films starts when characteristic FM-SC 
interaction energy e m exceeds the energy eo necessary for 
creation of a single SV. We have proposed experiments 
to detect the bound state of the magnetic dot and vor- 
tex and coupled pairs of magnetic and superconducting 
vortices. 
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